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Abstract. Periodic gravity-capillary waves propagating at a constant velocity at the surface of a fluid of infinite 
depth are considered. The surface tension is assumed to vary along the free surface. A numerical procedure is 
presented to solve the problem with an arbitrary distribution of surface tension on the free surface. It is found that 
there are many different families of solutions. These solutions generalize the classical theory of gravity-capillary 
waves with constant surface tension. An asymptotic solution is presented for a particular distribution of variable 
surface tension. 
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1. Introduction 

Gravity-capillary waves have very interesting properties. In particular there arc many different 
families of solutions characterized by dimples on their free surface. It is usual to study these 
waves by assuming that the surface tension is constant. The purpose of this paper is to show 
that the properties of gravity capillary waves remain qualitatively similar when the surface 
tension is variable along the free surface. 

Gravity-capillary waves with constant surface tension have been studied since the be¬ 
ginning of the century. The evidence of multiple solutions was first shown by Harrisson 
[1] and Wilton [2] who included surface tension in Stokes’s classical expansion for pure 
gravity waves. This work was later extended by Pierson and Fife [3], Nayfeh [4], Chen and 
Saffman [5] and others. Fully nonlinear numerical solutions were obtained by Schwartz and 
Vanden-Broeck [6], Chen and Saffman [7] and Hogan [8]. 

There arc physical situations in which the surface tension is variable along a free surface. 
For example the variation in surface tension might be due to a non-uniformity of the tem¬ 
perature or the presence of surfactants. Here we do not study how a particular distribution of 
variable surface might occur. We look at the general problem of periodic waves with variable 
surface tension. Then the equations force the surface tension to be a periodic function which 
can therefore be represented as a Fourier series. We derive a numerical procedure which 
incorporates this Fourier representation. In this approach all the coefficients in this Fourier 
series are free parameters whose values we can choose. We present numerical results for 
some typical values of these parameters. The results include as a particular case the findings 
of Vanden-Broeck [9] for pure capillary waves. 

Vanden-Broeck [9] showed that there arc symmetric and nonsymmetric waves with vari¬ 
able surface tension. For simplicity we restrict our attention to symmetric waves. We show that 
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gravity-capillary waves with variable surface tension have properties very similar to those with 
constant surface tension. Firstly, the linear theory predicts that there is a minimum speed below 
which periodic waves fail to exist. Secondly, there arc many different families of nonlinear 
waves. We show that this multiplicity of solutions arises, like in the theory with constant 
surface tension, from the fact that two linear waves whose wavenumbers arc multiples of each 
other can travel at the same speed. When the ratio of the wavenumbers is 2, the resulting 
nonlinear waves arc referred to as Wilton’s ripples. We present a semi-analytical theory of 
Wilton’s ripples with arbitrary variable surface tension. In addition we derive an asymptotic 
expansion for a particular distribution of surface tension. 

The formulation and the numerical procedure arc similar to those presented in Vanden- 
Broeck [9], A brief description is given in Section 2 and the reader is referred to Vanden- 
Broeck [9] for further details. The numerical results arc presented in Sections 3 and 4. The 
asymptotic expansion is derived in Section 5. An important question is how the gradient of 
surface tension is maintained. This is addressed at the end of Section 2, where we discussed 
the relation between the inviscid flows presented in this paper and the corresponding slightly 
viscous flows. 

2, Formulation and numerical procedure 

We consider a train of periodic waves with wavelength X traveling to the left at a constant 
phase velocity c. The fluid is inviscid, incompressible and of infinite depth. The flow is 
assumed to be irrotational. We introduce a frame of reference moving with the wave, with 
the origin at a trough and the y axis directed vertically upwards. As v —»■ — oo, the velocity 
approaches c. Gravity is acting in the negative y direction (see Figure 1). 

We introduce the complex potential / = 4> + i\|/ and choose / = 0 at x = 0. As discussed 
in the introduction, there are non-symmetric solutions. For simplicity we shall restrict our 
attention to waves which are symmetric with respect to x = y = 0. 

We denote by u and i; the horizontal and vertical components of the velocity. The basic 
equations can be written as (see for example Crapper [10] for a derivation) 


§xx + c l ) y>’ — 0, 

(1) 

<l>y = ° n >’ = hC*), 

(2) 

k(4>x 2 + <\>y 2 ) +gy+ —K = B on y = ri(x), 

(3) 

P 


4> —»■ cx as y —^ —oo, 

(4) 

4>(x + X, y) = cX + 4>(x, y), r\(x + X) = q(x). 

(5) 


Here g is the acceleration of gravity, T the surface tension, p the density, K the curvature 
and B the Bernoulli constant. The function y = x\(x) describes the shape of the free surface. 
Equation (1) is Laplace’s equation to be satisfied in the flow domain. Equations (2) and (3) 
are the kinematic and dynamic boundary conditions on the free surface. The conditions (5) 
impose the periodicity of the flow and (4) forces the flow to approach a uniform stream with 
constant velocity c as v —> — oo. 
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c 

-> 

Figure 1. Sketch of the flow and the coordinates. Only one wavelength of the wave is shown. As y —> —oo, the 
flow is uniform with constant velocity c. 


We assume a variable surface tension along the free surface. Equations (3) and (5) imply 
that T is then a periodic function of x, i.e. 

T(x + \) = T(x). (6) 


The problem defined by (1-5) is difficult to solve for two reasons. The first is that the 
boundary condition (3) is nonlinear. The second is that the shape of the boundaries (i.e. the 
function y = r](x)) has to be found as part of the solution. As shown by Stokes, this second 
difficulty can be avoided by taking ()> and \|; as independent variables. The flow domain is then 
mapped from the (x, y)-planc of Figure 1 to the lower half plane of the (c|), t|/)-planc. All the 
quantities x, y, u and v arc considered as functions of /. Relations (5) imply that there arc 
periodic in 4> with period c'k. From (6), it then follows that T is a periodic function of c|) with 
period c'k. The assumed symmetry of the wave about <)) = r|; = 0 implies that T is an even 
function of (j). We therefore represent T as a Fourier cosine series, i.e. 


T = 



(V) 


Here k = 2tt/X is the wavenumber and T 0 is a constant. 

This completes the formulation of the problem. We seek u — in as an analytic function of 
/, periodic in 4> with period c'k. This function must approach c as t|; —»■ —oo and satisfy (3) 
on r|/ = 0. 

To solve the problem numerically, we introduce dimensionless variables by choosing k as 
the unit length and c as the unit velocity. Next we define x — i0 by 

u — iv = e t-lH (8) 


Following Vanden-Broeck [9], we rewrite (3) as 


a 2t 


+ yy - a 


OO 

1 + a n cos(27T«(|)) 

n= 1 



where 



To 

pkc 2 


(9) 


and 


a = 


( 10 ) 
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Using the periodicity of the wave and the condition u — in —> 1 as i|r —— oo, we write 

OO 

x-i 0 = ^n„e- 2iM/ , (11) 

n =1 


where the coefficients u n arc real because of the assumed symmetry with respect to c|) = 0. 
Since u — iv = df/dz. 


— T+iO 


dx dy 1 

-b i— =-= e 

34> 3c|) u — iv 

Integrating (12) we rewrite (9) as 

a 2x 


e 2t M 

T + y / 0 


e T sin 0 dcj) — a 


1 + a n cos(2jtn(|)) 


n= 1 


t 30 „ 

e — = B. 

34> 


( 12 ) 


(13) 


Following Vanden-Broeck [9], we solve the problem numerically by truncating the infinite 
series in (13) after N — 1 terms. We find N + 2 unknowns u i, •••. m,v-i , Y- a an d B by 
collocation. We introduce the N mesh points 

^/ = ^7(/-l), 1 = 1,..., N. (14) 


By satisfying (13) at the mesh points (14), we obtain N nonlinear algebraic equations. 
The integral in (13) is evaluated numerically. An extra equation is obtained by fixing the 
steepness s of the wave (i.e. the difference of height between a crest and a trough divided by 
the wavelength). 

Using (12), we write this equation as 

e -T sin9d(|) = s, (15) 

where s is given. The integral in (15) is also evaluated numerically. 

To derive the last equation we first introduce the capillary number 


k = 4tt 2 


To 

PgT 2 


and the gravity parameter 


(16) 


|x = 2 tt^. (17) 

The parameters k, a and y arc related by the identity 


,a 

k = 4tt 2 -. 
Y 


(18) 


The last equation is given by (18) where the left-hand side is given. 

This system of N + 2 nonlinear equations with N + 2 unknowns is solved by Newton’s 
method. 

After finding the coefficients u n , the free-surface profile is obtained in the parametric form 
x = jc(4>), y = y(4>) by integrating numerically (12). We note that (7) defines T as a function 



The influence of variable surface tension on capillary-gravity waves 273 


of c|). This is an inverse formulation and T as a function of x is defined at the end of the 
calculations by (7) and x = jc (4>). 

For y = 0, the scheme reduces to the one used by Vanden-Broeck [9]. 

Let us mention that there is a variation of the numerical procedure in which the integral in 
(13) is avoided. The idea is to differentiate (13) with respect to c|). This leads 


e 2x — + ye T sin 6 —oi— 


3t 


3 


34 > 


34 > 


oo 

1 + a n cos(2jt«(|)) 

n= 1 



= o. 


(19) 


The numerical procedure is then the same but with this new equation instead of (13). Both 
numerical schemes were found to give equivalent results. 

Before concluding this section, we discuss how the inviscid solutions presented in this 
paper might approximate a slightly viscous flow (or more precisely a flow with high Reynolds 
number). For viscous flows, an extra boundary condition needs to be satisfied on the free 
surface, namely 


d T 

dx p'kc 2 


1 dv 3 u 

R ^ dx 3y ^ 


( 20 ) 


where 

c\ 

R = —. (21) 

v 

Flere v is the kinematic viscosity (see Lucassen-Reynders and Lucassen [11] for a derivation). 

For constant T and v f 0, (20) implies that the tangential stress vanishes on the free sur¬ 
face. For slightly viscous fluid and constant T, we can solve the flow problem approximately 
by first constructing an inviscid solution (not satisfying (20)). In a second stage, a boundary- 
layer solution is then constructed near the free surface where the tangential stress varies from 
its inviscid value 


1 / dv 3 u\ 
~R \ 3x + 3 y) 


( 22 ) 


to zero (see Batchelor [12], pp. 364). 

We propose to approximate the problem with variable surface tension in a similar way. 
First we solve the inviscid problem ignoring (20). These are the solutions presented in this 
paper. Then a boundary-layer solution can be added near the free surface where the tangential 
stress varies from its inviscid value (22) to the value 


dv fike 2 

predicted by (20). We note that our approach reduces to the classical one presented in Batche¬ 
lor [12] when T is constant. The boundary layer is not constructed here but the above argument 
suggests that the inviscid solutions calculated in this paper provide a good approximation for 
slightly viscous fluids away from a thin boundary layer near the free surface. 


3. Linear theory 

We can gain some preliminary insight into the problem by linearizing the equations around a 
uniform stream. 
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Assuming t 1 and 0 < 1 we linearize (19) as 

dx 0 [[ 1 301 

—- + Y0 - a— 1 + > a n cos(2tt/k|)) — =0. 

dcp d(b dtp 

_I n= 1 J _ 

When a n = 0 (i.e. for constant surface tension ), the problem has an exact solution 
x — i0 = Ae~ 2mf , p = 1 + k. 


( 24 ) 


(25a,b) 


Here A is a constant. 

Using (17) and (18), we can rewrite (25b) as 

c 2 _S} L+ T A 2 k 
2tt p X 

Equation (26) is the familial - dispersion relation of linear gravity-capillary waves. 
To describe the results further, it is convenient to introduce the new parameters 


X = 2jtK */ 2 = 


(26) 


(27) 


The parameters c and X can be viewed as the dimensionless velocity and the dimensionless 
wavelength if we choose (gT 0 / p) 1 / 4 as the unit velocity and (To/pg) ,/2 as the unit length. 
Using (27) we rewrite (26) as 


X 2tt 
2tt + X 


(28) 


Relation (28) is shown graphically in Figure 2 where we plot values of c 2 versus X (see 
solid curve). When X = 2 tt, c 2 has a minimum. 

When a n ^ 0 , we no longer have an exact solution because the linear problem has variable 
coefficients. Therefore we solve the linear problem numerically by the scheme of Section 2 
with (13) replaced by (24). The coefficients u n were found to decrease rapidly. For example 
u\ ~ 0-1047, u 2 ~ -0-0856, ..., « 7 g ~ 0-0000000174 for n g = 0-15, a/ — 0, j ^ 8 and 
k = 0-5. 

Numerical values of c 2 versus X are shown in Figure 2 for n g = 0-15 and a% = 0-45. As c/ g 
increases the curve moves down. Figure 2 shows that the relation between c and X for variable 
surface tension is similar to that with constant surface tension. In particular, the three curves 
of Figure 2 have a minimum. This behavior does not depend on the particular choice c/ g ^ 0, 
ci j = 0, /' ^ 8 and similar results were found for other values of a n . 

A consequence of the existence of minima in Figure 2 is that waves of different wave¬ 
lengths can travel at the same speed. In particular a wave of wavelength X and a wave whose 
wavelength is X/m (where m is an integer) can travel at the same speed if k = (2tt/X) 2 is 
equal to some critical value k„, . 

When the surface tension is constant (i.e. a„ = 0 for all n ), K m can be calculated analyti¬ 
cally as follows. Using (28), we see that waves of wavelength X and X/m travel at the same 
speed if 



(29) 
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Figure 2. Values of c 2 versus X. Here c and X are the di- Figure 3. Free surface profile of a Wilton ripple with 
mensionless velocity and the dimensionless wavelength D = A in (31). Here A = 0.1. Only one wavelength is 
defined in (27). The curves from top to bottom corre- shown, 
spond to ag = 0, ag = 0-15 and ag = 0-45 with ay = 0, 
j A 8. The three curves have a minimum. 

Relation (27) and (29) give 

Km = — - (30) 

m 

When the surface tension is variable (i.e. a n f 0), k„, cannot generally be calculated 
analytically. Therefore we present at the end of this section a numerical procedure to evaluate 
K m • 

For k = K m and constant surface tension, the solution (25 a) is incomplete and the general 
solution of the linear problem is 

x — i9 = Ae~ 2inf + De~ 2imnf , (31) 

where A and D are constants. The value of D cannot be determined within the framework of 
the linear theory. To find D, we need to use a nonlinear theory, i.e. use (19) instead of (24). 
For constant surface tension, Wilton [2] found D = ±A for k = k 2 = His calculation was 
based on an expansion in powers of the amplitude of the wave. There are two ‘Wilton ripples’ 
corresponding to D = A and D = — A. One has a crest dimple and the other a trough dimple. 
There are shown in Figures 3 and 4 for A = 0-1. 

For k = k,„ and variable surface tension, the general solution of the linear problem can be 
written as 

x —i6 = AL l (f) + DL 2 (f), (32) 

where L\(f) and L 2 (f ) are, respectively, the two linear solutions of wavelengths X. and \/m. 
We note that for constant surface tension L\(f) = and L 2 (f) = c~ 2l ™4, so that (32) 

reduces to (31). 

We conclude this section by describing the procedure to compute K m when the surface 
tension is variable . For simplicity we assume m = 2. We consider two linear waves whose 
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Figure 4. Same as Figure 3 with D = —A. Figure 5. Values of the critical capillary number k 2 ver¬ 

sus ag with a j = 0, j 7 ^ 8 . The variable surface tension 
is defined by (7). 


ratio of the wavelengths is equal to 2 (i.e. L \ (/) and L 2 (/)). In dimensionless variables they 
arc represented by the expansions (11) and 


■» = E 


/ -2i nnf 


u„e 


n =1 


(33) 


We require (11) to satisfy (24) and (33) to satisfy a modification of (24) (denoted by (24')) 
obtained by replacing with y and a by the new parameters y and a. These new parameters 
satisfy 

a a 

— = 4— and ya = ya. (34, 35) 

Y Y 

Relations (34) and (35) express the facts that the ratio of the wavelengths is equal to 2 and 
that the waves travel at the same speed. We truncate the infinite series in (11) and (33) after 
N terms and satisfy (24) and (24') at the mesh points (14). This yield 2N equations for the 
2N + 4 unknowns u i, n 2 > ..., un, u[, u' 2 , ..., u' N , a, y, a and y. Relations (34) and (35) 
provide two more equations. We obtained the last two equations by fixing the steepness of 
the waves. We note that the numerical results are independent of the actual values chosen 
for the steepness since the problem is linear. This system of algebraic equations is solved by 
Newton’s method. The value of k 2 is then given by 


4jt 2 o( 

K2 = - 

Y 


(36) 


Values of k 2 versus a s with a j = 0, / ^ 8 arc shown in Figure 5. 

We note that the scheme gives in addition to the value of k 2 , the two linear solutions L \ (/) 
and L 2 (f ). These solutions will be useful for the theory of Wilton’s ripples presented in the 
next section. 
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Figure 6. Free surface profiles for ag = 0-15, aj = 0, j ^ 8 and s- = 0-0032 (Figure 6a), .y = 0-0153 (Figure 6b) 
and s = 0-028 (Figure 6c). Flere s is equal to the ordinates of the crests. The solid curves are the nonlinear 
solutions and the broken curves correspond to the right-hand side of (37). 



x 

Figure 6. Continued. 
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Figure 7. Free surface profiles for k = 0-33, cij — 0, 
) / 8 and ag = 0 (solid curve), ag = 0-15 (broken 
curve) and ag = 0-30 (dotted line). 



4. Nonlinear theory 

We used the numerical scheme of Section 2 to compute solutions for various values of s, 
a n and ic. The coefficients u n were found to decrease rapidly as n increases. For example 
~ 0-02239, ui ~ -0-0162, ..., m 78 *=» 0-00000002 for o 8 = 0-15, aj = 0, j ^ 8, k = 0-5 
and 5 = 0-007. 

First we checked the scheme by reproducing the results of Schwartz and Vanden-Broeck 
[6] for constant surface tension . Next we used the scheme to generalize the Wilton ripples 
(i.e. waves for k = k 2 ) to the case of variable surface tension. Here we present typical results 
for «g / 0 and a j = 0 for j f 8. For a given values of a 8 , the corresponding value of k 2 can 
be calculated by the scheme described at the end of Section 3 (see Figure 5). We then compute 
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a nonlinear solution for these values of c/ x and k by the scheme of Section 2. We denote the 
computed values of t — i0 by Nl(f ). In the limit as s —> 0, we expect that 

ALdf) + DL 2 (f), (37) 

where L\{f) and L 2 (f ) arc the two Iinear solutions calculated at the end of Section 3. To 
check (37) we find A and D by imposing two conditions. The first is that Nl(f ) and AL\ (/)+ 
DL 2 (f) have the same steepness. The second is that 

X/2 

[2V/ (4>) — ALi(fy) — B L 2 (<\>)] 2 d§ (38) 

is minimum. The last condition imposes that the difference between the right- and left-hand 
sides of (37) is minimum in the L 2 norm. We can now compare the right- and left-hand sides 
of (37). Typical results for = 0-15, cij = 0, j / 8 arc shown in Figure 6. The solid curve 
corresponds to Nl(f ) and the broken curve to AL\( f ) + DL 2 (f). The results show that the 
agreement between the solid and broken curve improves as the steepness s is decreased. Since 
we chose the origin at the through, the steepness is simply the ordinate of the crests. This 
constitutes a check on (37) and on the various numerical procedures. 

Finally we present in Figure 7, numerical results for k = 0-33 and various values of a g , 
ctj = 0, j i=- 8. Here we found that the dimples arc less pronounced and that the wave is 
getting flatter near - the trough as we increase a%. 


5. Perturbation results 


In the previous sections, the surface tension is a prescribed (periodic) function of (j). The 
numerical procedures can be extended to cases in which the surface tension is prescribed as a 
function of other quantities along the free surface. As an example we present in this section 
results for a variable surface tension which depends on the curvature of the free surface. More 
precisely we replace (7) by 

T = Toll - m, (39) 


where e is a given constant. The numerical procedure of Section 2 applies unchanged if we 
replace the square bracket in (9) by 


, e 30 
1 + — e — 
k 34> 


(40) 


We denote this equation by (9') The linearization of (9') is (24) with a„ = 0 for all n. In other 
words the linearization of (9') is the classical linear - problem with constant surface tension 
whose solution is given by (25 a) and (25b). This nice property suggests that for the variable 
surface tension (39), it should be easy to improve on the linear - theory by constructing an 
expansion in powers of the amplitude of the wave. This is achieved up to second order in this 
section. This provides us with some analytical results to supplement the numerical calculations 
presented in the previous sections. 

Since we are only calculating solutions up to second order in the amplitude, it is more 
convenient to use Cartesian coordinates instead of the variables 4> and r|r. 
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The governing equations in dimensional variables arc then (1), (2), (4), (5) and 


1 9 9 2o 

2 (^ + ^ v ) + gy — 


p L 


1+6 




(1+T1.V 2 ) 372 . 


4a 


1 


d + v+ = + on y = nW ' <41) 


Here we chose the origin of y such that the Bernoulli constant is c 2 /2. 

We seek a solution in the form of an expansion in powers of a small parameter e. Thus we 
write 


riCO = 6ru(x) + 6 2 r) 2 <» 4-, 

4>(x, y) = cx + €<(>i(. x , y ) + (r§ 2 (x, y) + 

c = Co + 6Cl + e“C 2 + • • •, 


(42) 

(43) 

(44) 


The parameter 6 is a measure of the wave amplitude. Here we define it as e = ak where a is 
the first Fourier coefficient of q(x). 

Substituting (42)—(44) in (1), (2), (4), (5) and (41) and equating the coefficients of e, we 
obtain 


4+M + <t>lyy — o, (45) 

hy = com* on >’ = 0, (46) 

T 

co<k* + gm-m.« = o, (47) 

P 

4>i —^ 0 as y —> — oo, (48) 

<t>t(* + K y) - <t>i(^» y), m(* + X) = mC*)- (49) 

Similarly, equating the coefficients of e 2 , we obtain 

4>2** + §2 yy =0- (50) 

<t> 2 y-C 0 ri 2 ^ =-<|)iyyTii + (Cl+ c|)u)riu on y = 0, (51) 


r ° ^ X ^ 

- hlxx + < 4 ) 92 * + gT\2 — 

P 

T 

-Co^ixyhi - Ci$u - j((|m T 2 + 6t>iy 2 ) + —6T] 2 i.„ on y = 0, 

p 

4> 2 —^ 0 as y — ► — oo, 

<te(* + k, y) - <\> 2 (x, y ), r\ 2 (x + X) = r\ 2 (x). 

We note that (45-49) is the I i near problem with constant surface tension. 

We first assume that 

1 

k # — • 

m 


(52) 

(53) 

(54) 


( 55 ) 
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As shown in Section 3, the solution of the linear problem (45—49) is unique and given by 


<|)i(x, y) = — coA sin ( kx ), 


(56) 

T| 1 (x) =AiCOS(fcx), 


(57) 

2 § , r 0 , 

Co — 7 H- 

p 


(58) 

where k = 2tt/X. is the wavenumber and Aj is 
the wave. The definition of e implies that A i = 
Next we use separation of variables to write 

a constant. Here we chose x = 0 at a crest of 
1 /k. 

the solution of (50), (53) and (54) as 

OO 

(b(x, y) = ^ F n e l,ky sin(nkx) 

n= 1 


(59) 

and we seek r\ 2 (x) in the form 



OO 

riaCO = Eq + ^ E n cos(nkx). 


(60) 


n =2 


We started the sum in (60) at n = 2 in accordance with our definition of e 

Substituting (59) and (60) in (51) and (52) and equating the constant terms in (51) and the 
coefficients of sinx and cosx in (51) and (52), we obtain 

T 

Eq = ——eA^k 4 , kFi = —kc\A\, kc^Fx = kcoC\Ax. (61,62,63) 

2 pg 

Similarly, equating the coefficients of sin(2x) and cos(2x) in (51) and (52) yields 

2 kF 2 + 2kc 0 E 2 = k 2 c 0 Ax 2 , (64) 

4 k 2 F 2 — + gE 2 + 2 kc 0 F 2 = ifcVV + l -\k 4 Ax 2 . (65) 

P - 2 p 

Finally, by equating coefficients of sin(nx) and cos(nx) in (51) and (52) for n — 3, 4, • • • 
we have 


nkF n + nkcoE n = 0, 


( 66 ) 


2 Tq 

n 2 k~E n -1- gE n + nkcoF n = 0. 

P 

Solution of (62-65) gives 
F\ = 0, ci = 0 
and 

& = 1 c ° 2 ~ jfg p = C « T «, 3 ~W 1 
- 2t(2ta-f)' 2 P2 *a-{- 


(67) 


( 68 ) 


(69, 70) 
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Here we have used A\ = 1 /k and (58). 

Eliminating F n between (66) and (67), we get 


( nk) 2 -1- g — nkcf 

P 


E n = 0. 


Relations (66) and (67) imply 

E n — F n — 0 for n — 3,4,..., 


(71) 


(72) 


provided the square bracket in (71) is different from zero. This requirement is satisfied if (55) 
holds. 

Up to second order the solution is defined by (42-44), where the various terms arc given by 
(56-61), (68-70) and (72). We note that the denominators of the second-order approximation 
(69), (70) vanish when k = Similar singularities in the higher-order terms occur when 
k = I /m. These singularities arc related to the fact that the solution (56), (57) of the linear 
problem is incomplete when \c — l/m. 

We now examine in details the case k = 1/2 (i.e. we develop a theory of the Wilton ripples 
for the variable distribution of surface tension (39)). 

The complete solution of the linear - problem is then 

4>i (jc, y) = — c 0 Aie^sin (kx) — coA 2 e 2<:> ’sin (2 kx), (73) 


T]i(x) = Aicos (kx) + A 2 cos (2 kx), (74) 

where Co is defined by (58). Relations (73), (74) are equivalent to (31). The difference is that 
the independent variables are 4> and r|/ in (31) and x and y in (73), (74). Following the analysis 
in the case k f I//;/, we substitute (73) and (74) in the right-hand sides of (51) and (52) and 
we seek 4> 2 and r\ 2 in the form (59) and (60). Equating the constant terms and the coefficients 

of sinx, cosx, sin2x, cos2x.in (51) and (52) give equations for A 2 , c i, F„ and E n . Here 

we present only the calculations leading to the determination of the constants A 2 and c\ in the 
linear - solution. 

By equating the coefficients of sin x in (51) and of cos x in (52), we have 


kF\ = jk 2 coAiA 2 — kc\A\, (75) 

T 

kcoFi = \k 2 C(f A\A-> + kc$c\ A\ -f- 4— ck^A^Ai. (76) 

P 

Similarly, by equating the coefficients of cos 2x in (52) and sin2x in (51) we have 

2kF 2 + 2 kc Q E 2 = k 2 c Q Ai 2 - 2 kc { A 2 , (77) 

T 1 T 

4k~E-> -b gE-> T 2kco F 2 = — k" cq“ A\" + 2kcoC\A-> -f- i— ck^Ap. (78) 

P 2 P 

Solving (75-78) we get 

i, A r.To~ , 3 A 2 
ci = IjkcoAo — 2 —ck — 

P c 0 


( 79 ) 
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and 


A 


2 

2 







Q) 


To „ 9 

4— Ik 2 


( 80 ) 


Relation (80) shows that Ac, increases as e increases. Similarly, Equation (79) implies that |ci | 
decreases as e increases. 

Note that, when e = 0, we recover the relation 


A 2 = ±\A r (81) 

obtained by Wilton [2] for constant surface tension. The corresponding free-surface profiles 
arc shown in Figures 3 and 4. 


6. Conclusions 

We have generalized the theory of periodic gravity capillary waves by allowing the surface 
tension to vary along the free surface. We have found that the linear theory predicts a minimum 
of the phase velocity. We have shown numerically and analytically that this minimum implies 
the existence of multiple solutions similar to the ones obtained for constant surface tension. 
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